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ABSTRACT 

We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the 
M 81 group. The far-infrared data have been obtained as part of the 'Spitzer Infrared Nearby Galaxies 
Survey SINGS\ Maps of the distribution of atomic hydrogen (HI) have been obtained through 'The 
Hi Nearby Galaxy Survey THINGS'. The Spitzer observations provide a first glimpse of the nature 
of the non-atomic ISM in these metal-poor (Z~0.1Zq), quiescent (SFR~0. 001-0. 1 Mq yr~^) dwarf 
galaxies. Dust emission is detected in five out of the seven targets (the two systems with the lowest 
star formation rates are non-detections). Most detected dust emission is restricted to Hi column 
densities > 1 x 10^^ cm~^ and almost all regions of high Hi column density (> 2.5 x lO^'^cm"^) 
have associated dust emission. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 
and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic 
features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral 
galaxy by a factor ot '--^ 7 . The aromatic features in Holmberg II (which has only a slightly lower 
gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises 
that the strength of the aromatic features is not a simple linear function of metallicity. Whereas the 
H I masses are well-constrained, model dependencies make it difficult to measure the dust masses with 
a high degree of confidence. We estimate dust masses of ^10"'-10® MQfor the M81 dwarf galaxies, 
resulting in an average dust-to-gas ratio (Mdust/Mni) of '--^ 3 x 10^^ (1.5 x 10^"^ if only the Hi that 
is associated with dust emission is considered); this is an order of magnitude lower than the typical 
value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation 
rate at 70/im as compared to the more massive galaxies in SINGS by a factor of ^ 2. However, the 
average 70/im/160/xm ratio in the sample dwarf galaxies is higher than what is found in the other 
galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in 
conjunction with a higher dust temperature in the dwarfs (likely due to the harder radiation fields in 
the low metallicity environments). 

Subject headings: galaxies: dwarf — galaxies: irregular — galaxies: ISM — infrared: galaxies — 
galaxies: individual (Ho II, M 81 dwA, DDO 53, Ho I, M 81 dwB, IC 2574, DDO 165) 
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1. INTRODUCTION 

Nearby dwarf galaxies have proven to be ideal labora- 
tories to investigate how stars form out of gas and how, 
in turn, violent star formation shapes the ambient inter- 
stellar medium (ISM). These systems are highly suscep- 
tible to the formation of shells and holes in the neutral 
gas phase; more energetic star formation can lead to the 
formation of gaseous outflows from these systems. Previ- 
ous studies have shown that dwarf galaxies can be used 
as testbeds of 'simple' prescriptions for star formation 
(since these systems are typically in solid-body rotation 
and are therefore less affected by shear in the ISM) and 
for understanding the connection between mechanical en- 
ergy input into the ISM ('feedback') and future star for- 
mation. Given their low metallicities, they also provide 
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a unique opportunity to study the eonditions of the ISM 
in environments that may resemble those in the earUest 
starforming systems at high redshift. 

Based on their ISM properties, dwarf galaxies can be 
roughly divided into two subgroups: gas-rich dwarf irreg- 
ulars (dirrs) and gas-poor dwarf spheroidals/ellipticals 
(dSph/dE). The gas-rich dIrrs are particularly interest- 
ing for studies of current star formation, since they still 
contain the 'fuel' for star formation. Though many stud- 
ies have shown that these objects are rich in atomic 
hydrogen (Hi; for a review see Skillman 1996), little 
is known about their molecular gas properties. Many 
searches for molecular gas (through observations of the 
most abundant tracer molecule, CO) have been per- 
formed but very few dwarf galaxies have been detected in 
CO so far: no dwarf has been detected in CO at metal- 
hcities < 10% Zq (Taylor et al. 1998; Barone et al. 2000; 
Leroy etal. 2005). 

Similarly, little is known about the dust properties in 
faint, low metallicity dwarf galaxies. Some of the bright- 
est dwarfs have been detected with previous far-infrared 
(far-IR) observatories such as IRAS and ISO (see, for ex- 
ample. Hunter etal. 1989, Gallagher etal. 1991, Melisse 
& Israel 1994a,b, Hunter etal. 2001). These early stud- 
ies found that dwarfs typically have higher dust tem- 
peratures than those derived for more massive galaxies, 
with the peak of the infrared spectral energy distribution 
(SED) shifted to shorter wavelengths. 

The more typical 'quiescent' dirr galaxies (with 
SFR<0.1 Mq yr~^) had remained undetected in the far- 
IR prior to the advent of the Spitzer Space Telescope; 
its dramatically increased sensitivity compared to pre- 
vious observatories has opened up the low-metallicity 
regime of the extragalactic ISM to exploration in the 
far-IR. Early Spitzer observations have naturally concen- 
trated on the brightest and most extremely metal-poor 
dwarfs. Houck etal. (2004) show that SBS 0335-052, 
one of the most metal-poor dwarf galaxies known, has 
an exceptional SED that is shifted blueward with a peak 
at ~ 28fim (compared to ~80 /xm for more metal-rich 
star burst galaxies and >100 /um for local star-forming 
spirals). Engelbracht etal. (2004) presented observa- 
tions of the dwarf galaxy NGC 55, and a larger sam- 
ple of galaxies showed that there is a metallicity thresh- 
old above which emission from aromatic features appear 
(Engelbracht etal. 2005). Smith et al. (2006) used 
Spitzer spectroscopy to study how the spectral signatures 
of aromatic features change as a function of metallic- 
ity. Other recent Spitzer studies of actively star forming 
dwarf galaxies include spectroscopic observations by Wu 
et al. 2006 and O'Halloren et al. 2006 and imaging studies 
of star-forming dwarfs and Local Group galaxies (includ- 
ing dwarf galaxies) by Rosenberg et al. 2006 and Jackson 
et al. 2006. Cannon etal. 2005, 2006a, b have presented 
detailed Spitzer studies of some of the most luminous 
and nearby dwarfs (IC 2574, NGC 1705, and NGC 6822). 
These case studies have shown that localized star forma- 
tion has dramatic effects on the multiphase ISM, altering 
the relative strengths of nebular, far-IR dust, and radio 
continuum emission. 

Here we discuss the distribution of dust and atomic 
hydrogen in seven dwarf irregular galaxies in the M81 
group of galaxies. Our sample spans two orders of mag- 
nitude in Hi masses and a similar range in star forma- 



tion rates, from < 0.001 M0 yr^^ (Ha non-detections) 
to ~O.1M0 yr"^. The general properties of the sample 
dwarfs are summarized in Table 1. All these dwarf ir- 
regular galaxies are part of the ' Spitzer Infrared Nearby 
Galaxies Survey' {SINGS, see Kennicutt et al. 2003) and 
have also been included in 'The Hi Nearby Galaxy Sur- 
vey' {THINGS, Walter etal. 2005). 

This paper is organized as follows: in § 2 we summarize 
the Spitzer MIPS and the VLA Hi observations; in § 3 
we present our results, i.e. the distribution of the atomic 
gas and the dust, IRS spectra of two galaxies, estimates 
of the dust masses and a comparison to other galaxies in 
the SINGS sample. In § 4 we present a summary of our 
study. 

2. OBSERVATIONS 
2.1. Spitzer MIPS Observations 

The 24, 70, and 160 /um data were taken with MIPS 
on the Spitzer Space Telescope as part of the SINGS 
survey (Kennicutt et al. 2003). The observations were 
obtained using the scan-mapping mode in two separate 
visits to each galaxy (facilitating removal of asteroids and 
detector artifacts). Each pixel in the map was observed 
40, 20, and 4 times at 24, 70, and 160y^m, respectively, 
resulting in integration times per pixel of 160, 80, and 
16 s, respectively. All MIPS data were processed using 
the MIPS Instrument Team Data Analysis Tool (Gordon 
etal. 2005). Systematic uncertainties (e.g., detector non- 
linearities, time-dependent responsivity variations, back- 
ground removal, etc.) limit the absolute flux calibration 
to ~4%, 7% and 12% in the MIPS 24iim, 70 ^m and 
160 Aim bands. The FWHM of the MIPS PSFs are 6", 
18", and 40" at 24/zm, 70/zm and 160 /xm, respectively. 
For more details on the MIPS data reduction see Bendo 
et al. (2006) - for a general description of the SINGS 
observing strategies see Kennicutt et al. (2003). 

The flux densities presented in this paper were de- 
rived for apertures much larger than the MIPS PSFs, 
and therefore aperture corrections have not been ap- 
plied. The M81 group is located in a direction where the 
Galaxy is rich in infrared cirrus (de Vries et al. 1987). In- 
spection of individual images reveals that Galactic cirrus 
emission is indeed present in the longer wavelength MIPS 
images. However, this emission is distributed over much 
larger angular scales than the sources of interest and can 
thus be easily separated from the galaxies presented in 
this study. 

The aperture for each individual galaxy has been cho- 
sen carefully to encompass all the emission visible in all 
three MIPS bands; apertures were compared with IRAC 
band 1 and Hi imaging to ensure that the total galaxy 
extent (i.e., gas and stars) was measured. For each indi- 
vidual galaxy, we used the same aperture to extract the 
flux densities from the three MIPS bands. To account 
for the variations in the background (both instrumen- 
tal and due to Galactic Cirrus) we have defined mul- 
tiple background regions for each galaxy that contain 
the same area as the target aperture. The (background 
subtracted) source flux densities we derived using this 
technique are the same (within the errors) as the values 
derived by Dale et al. 2005, 2006 in their study of the en- 
tire SINGS sample. For ccmsistcmcy we therefore adopt 
the values of Dale et al. 2006 for our study. The global 
flux densities are summarized in Table 1 — note that 
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TABLE 1 

Properties of the sample dwarf galaxies in the M 81 Group 



Galaxy 


j^a 

(Mpc) 


Fhi 
(Jykms~^) 


Mm 
(108 Mo) 


log(F(Ha))'' 
(ergs~^ cm~^) 


12+log(0/H)<: 


S(24/im)d 

(Jy) 


S(70/im)<i 

(Jy) 


S(160Aim)'^ 

(Jy) 


IC 2574 


4.02 


386.7 


14.75 


-11.27 


7.94±0.06 


0.28±0.013 


5.55±0.42 


11.75±1.50 


Holmberg II 


3.39 


219.3 


5.95 


-11.30 


7.68±0.03 


0.20±0.008 


3.67±0.26 


4.46±0.58 


Holmberg I 


3.84 


40.1 


1.40 


-12.43 


7.54±0.34 


0.013±0.002 


0.41±0.08 


0.90±0.17 


DDO 165 


4.57 


35.0 


1.72 


-12.93 


7.76^10.18^* 


<0.014 


<0.15 


<0.33 


DDO053 


3.56 


20.0 


0.60 


-12.24 


7.77±0.1 


0.029±0.001 


0.40±0.03 


0.50±0.11 


M 81 DwB 


5.3 


3.8 


0.25 


-12.82 


7.85±0.17 


0.009±0.001 


0.15±0.03 


0.39±0.18 


M 81 DwA 


3.55 


4.1 


0.12 






<0.002 


<0.17 


<0.15 


^Distances arc 
^Ha fluxes arc 
Characteristic 


from Karachentscv et al. 2002, 2003 
from Kennicutt et al. 2006, Lcc 2006 
gas phase oxygen abundances from 


Moustakas et al. 2007 


based on the 


Pilyugin & Thuan 


2005 strong 


-line calibration of 



R23=([OII] + [OIII])/H/3; 12+log(O/H)0=8.7 (Asplund et al. 2004). 
^Flux densities adopted from Dale et al. 2006. 

^Based on the luminosity— metallicity relation, Moustakas et al. 2007 




Fig. 1. — IC 2574: (a) integrated THINGS Hi map (contour shown at Nhi=10^^ cm"^ in all panels); (b) R-band image; (c) IRAC 
3.6/im image; (d), (e), and (f): 24 ^m, 70 /im and 160 /im images, respectively. The sizes of the Hi and MIPS beams are given in the lower 
left corners of their respective panels. 



M 81 dwA and DDO 165 are MIPS non-detections. The 
reader is referred to Dale et al. 2006 for the IRAC flux 
densities of the dwarf galaxies in our sample. 

2.2. VLA Observations 

Hi data for six of the seven M81 group dwarfs pre- 
sented here were obtained as part of 'The Hi Nearby 
Galaxy Survey' (THINGS), a survey to obtain high- 
resolution NRAO^^ VLA H I imaging for 35 nearby galax- 
ies (Walter et al. 2005). H i data for DDO 165 were taken 
from Cannon ct al. (in prep.). 

For THINGS, each galaxy was observed with the VLA 
in D, C and B configurations with typical integration 
times of 1.5 hours, 2.5 hours and 7 hours, respectively. 
The calibration and data reduction were done using the 
AIPS package^^. The absolute flux scale for the data was 
determined by observing the quasar 3C286 in all observ- 
ing runs. The time variable phase and amplitude calibra- 

The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 

The Astronomical Image Processing System (AIPS) has been 
developed by the NRAO. 



tion were done using the nearby, secondary calibrators 
1313-1-549 and 1252-1-565 which are unresolved for the ar- 
rays used. The uthdata were inspected for each array and 
bad data points due to cither interference or cross-talk 
between antennae were removed, after which the data 
were calibrated. After flnal editing, all data for each tar- 
get were combined to form a single dataset which was 
subsequently used to create maps of the brightness dis- 
tribution on the sky as a function of frequency /velocity 
(data cubes). 

In order to remove the continuum from the line data 
we flrst determined the line-free channels in our obser- 
vations and subtracted the continuum emission in the 
Mi>-planc. After that, datacubes (1024 x 1024 pixels x 
80 channels each) were produced using the task IMAGR in 
AIPS. To boost the angular resolution while still main- 
taining a reasonable noise, we use a robust parameter 
of 0.5 for the final imaging. To ensure that we reach iden- 
tical beam sizes for all THINGS observations, the data 
were subsequently convolved to a common resolution of 
10" . This resulted in a typical rms noise per channel of 
0.5 mJybeam"^ for a 2.5 km s~^ channel (corresponding 
to NHi=1.5xl0^^cm~^). To separate real emission from 
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Fig. 2. — Holmbcrg II: (a) integrated THINGS Hi map (contour shown at Nhi=10^^ cm^'^ in all panels); (b) R— band image; (e) 3.6/im 
image; (d), (e), and (f): 24 //m, 70 /xm and 160 //m images, respectively. The sizes of the Hi and MIPS beams are given in the lower left 
corners of their respective panels. 
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Fig. 3. — Holmberg I: (a) integrated THINGS Hi map (contour shown at Nhi=10^^ cm"^ in all panels); (b) R-band image; (c) IRAC 
3.6/nm image; (d), (e), and (f): 24 /im, 70jUm and 160 /im images, respectively. The sizes of the Hi and MIPS beams are given in the lower 
left corners of their respective panels. 



noise in the final integrated H I maps, we only consider re- 
gions which show emission in consecutive channels above 
a set level 2(t) in slightly convolved (20") cubes. Note 
that the data for DDO 165 arc at a resolution of 21". 

The fluxes in the integrated THINGS H i map arc cor- 
rected for the fact that typically the residual flux of the 
source in cleaned channel maps is overestimated (some- 
times by a factor of a few) due to the different beam 
sizes of the dirty and cleaned beams (for details sec, e.g., 
Jorsater & van Moorsel 1995, Walter & Brinks 1999). In 
the integrated Hi maps this typically leads to flux cor- 
rections of order 25-40%. To correct for this, we have 
scaled the residual fluxes by the ratio of the dirty and 
clean beam sizes and estimate that our column densities 
are correct within 10% (including the intrinsic uncertain- 



ties of the flirx calibration; 
given in Table 1). 



the inferred Hi masses are 



3. DUST AND Hi CHARACTERISTICS 

In Figs. 1-7 we present images of the individual galax- 
ies (in order of decreasing Hi mass). For each galaxy, 
we show six panels: (a) is the integrated THINGS Hi 

map at 10" resolution (only exception: DDO 165 beam- 
size: 21"); unless otherwise stated, one Hi contour is 
drawn at Nhi=1 x 10^^ cm^^ (i.e., close to the canonical 
star formation threshold, e.g., Skillman 1996). An opti- 
cal R-band and the SpitzerlRAC band 1 (3.6 /xm) image 
of the galaxies are shown in panels (b) and (c) (the op- 
tical images have been observed either at the Calar Alto 
2.2 m telescope or are taken from the ancillary SINGS 
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Fig. 4. — DDO 165: (a) integrated Hi map at 21" resolution (data from Cannon et al. in prep.; contour shown at Nhi=10^^ cm^^ in all 
panels); (b) R— band image; (c) IRAC 3.6/im image; (d), (e), and (f); 24 /im, 70 /^m and 160 /im images, respectively. The sizes of the Hi 
and MIPS beams are given in the lower left corners of their respective panels. 
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Fig. 5. — DDO 53: (a) integrated THINGS Hi map (contour shown at Nhi=10^i cm~^ in all panels); (b) R-band image; (c) IRAC 
3.6//m image; (d), (e), and (f): 24 /^m, 70 ^im and 160 //m images, respectively. The sizes of the Hi and MIPS beams are given in the lower 
left corners of their respective panels. 



data archive). Panels (d), (e) and (f) are the MIPS 24, 
70 and 160 fim images of the same area. All panels show 
the same H i contour as presented in the first panel. The 
bcamsizes for both the H i and MIPS images are shown 
in the lower left of the respective images. 

For each galaxy with THINGS H i imaging, we com- 
pare the radial profiles of the Hi, 70/Ltm and 160 /xm 
MIPS images in Fig. 8 (as discussed for the individual 
systems below). The deprojection parameters (inclina- 
tion, position angle) used for the creation of the radial 
profiles have been derived from the H i maps (see caption 
Fig. 8). 

In the following, we briefly discuss the individual galax- 
ies: 

IC 2574 (Fig- 1): IC 2574 is the largest galaxy in our 
sample and its H I morphology is dominated by the pres- 
ence of Hi holes (Walter & Brinks 1999). The brightest 
region in the MIPS bands is the supergiant shell (SGS) 
region in the north-east (Walter etal. 1998); a spatially 



resolved Spitzer case study of this SGS region is presented 
in Cannon et al. 2005. The elevated emission toward 
the south-east corner in the 160 /xm image is caused by 
Galactic Cirrus emission (but this emission can be sep- 
arated from the emission of IC 2574). Dust emission 
traced by the 70 ^m emission is detected out to galacto- 
centric radii of 7' (~7kpc, see Fig. 8). 

Holmberg II (Fig. 2): The distribution of Hi in Holm- 
berg II is also characterized by the presence of numerous 
Hi holes (Puche etal. 1992). As in the case of IC 2574, 
the changing morphologies in the individual MIPS bands 
and the corresponding changing spectral energy distribu- 
tions stress the importance of local eff'ects in character- 
izing the far-infrared emission (e.g.. Cannon etal. 2005). 
The radial surface brightness profiles (Fig. 8) show that 
dust is detected out to at least 4' ('^4kpc). 

Holmberg I (Fig. 3): The H I distribution in Holmberg I 
is characterized by one giant Hi hole (Ott etal. 2001). 
The Hi structure encompasses the optical emission and 
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Fig. 6. — MSldwB: (a) integrated THINGS Ki map (contour shown at Nhi=10^^ cm"^ in all panels); (b) R-band image; (c) IRAC 
3.6/im image; (d), (e), and (f): 24 /im, 70 /um and 160 /im images, respectively. The sizes of the Hi and MIPS beams are given in the lower 
left corners of their respective panels. 
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Fig. 7. — MSldwA: (a) integrated THINGS Ul map (contour shown at Nhi=3x10^° cm"^ in all panels); (b) R-band image; (c) IRAC 
3.6/im image; (d), (e), and (f): 24fim, 70 /im and 160 /im images, respectively. The sizes of the Hi and MIPS beams are given in the lower 
left corners of their respective panels. 



faint star formation is present on the rim toward the 
south-east (see also Sec. 3.4). Tliis is the region where 
faint emission is detected in all three MIPS bands. At 
70 /zm, there is also diffuse emission present toward the 
western H i rim (though this emission is of very low S /N) . 

DDO 165 (Fig. 4): DDO 165 shows extended emission 
in both the H i and the optical, but is not detected in the 
MIPS bands. The brightest emission seen in the 70^m 
and 160/xm bands (toward the north) is a background 
galaxy (SDSS J130639.44+674456.4 at z=0.139). There 
is some emission seen at 24 iim and 70 /im toward the cen- 
tre of DDO 165 (not coincident with the peak of the H i 
emission, and outside the main optical body of DDO 165) 
- future, higher sensitivity observations (MIPS and H l) 
are needed to see if this emission is indeed physically re- 
lated to DDO 165. Note that DDO 165 has one of the 
lowest star formation rates in our sample. 

DDO 53 (Fig. 5): The Hi distribution shows two 
peaks, and the galaxy is detected in all three MIPS 
bands: The brightest emission seen at 24 fiva. is associated 



with the northern and southern H i peak; the northern re- 
gion is also the strongest in the 70 ^m image. The 160 jim 
data are noisy, but 160 yum emission is still present in the 
regions seen in the 70 /zm image. The compact nature of 
DDO 53 is also evidenced by the radial surface brightness 
profiles shown in Fig. 8. 

MSldwB (Fig. 6): MSldwB is the galaxy with the 
lowest measured star formation rate in our sample. The 
galaxy shows a compact structure in Hi and the MIPS 
bands (cf. Fig 8). The detection at 160 /im is marginal 
(sec Tab. 1) and is surrounded by elevated background 
emission present. 

M81 dwA (Fig. 1): MSldwA is the faintest dwarf in 
our sample. Similar to Holmberg I, the Hi distribution 
is characterized by one large H i shell which encompasses 
most of the optical galaxy. No ongoing star formation 
has been detected in MSidwA (Miller & Hodge 1994); 
it is also a MIPS non-detection. This may be explained 
by the fact that H i column densities do not reach values 
higher than Nhi=5x 10^*^ cm~^ in this galaxy (i.e. the 
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Fig. 8. — Comparison of the Hi (solid), 70 /xm (sliort daslied) 
and 160 /im (long dashed) radial profile of five M81 group dwarf 
irregular galaxies. The following deprojection parameters (posi- 
tion angle, PA; inclination i) where used to create these profiles: 
IC 2574 [PA=-|-60°, i=57°], Holmberg II [PA=-15°, i=26°]. Holm- 
berg I [PA=55°, i=46°], DD053 [PA=-45°, i=40°], MSldwB 
[PA=-60°, 1=46°]. 

Hi column densities are below the canonical threshold 
for star formation). 

3.1. HI Threshold for warm dust 

From an inspection of the morphologies (Figs. 1-7) 
and radial profiles (Fig. 8) it is clear that the dust emis- 
sion appears to be related to the distribution of the H i, 
at least to first order. To quantify this, we investigate 
if there is a certain Hi threshold above which most of 
the dust emission is present. In the following, we use 
the 70 /im data as a tracer for the warm dust emission 
as they are (unlike the 24 /im data) not affected by the 
presence of contaminating point sources (stars and back- 
ground objects) in the field. Furthermore, they have 
higher resolution and signal-to-noise than the 160 /im 
measurements. In Figure 9, we plot histograms of the 
distribution of the 70 /zm flux density above a threshold 
of 1.8 MJysr"^ (about the 5a level, i.e., encompassing 
most of the detected dust emission) as a function of H i 
column density for each galaxy (solid histogram). For 
each galaxy we also show the pixel-by-pixcl distribution 
of all H I column densities (dashed histogram) . From 
Fig. 9 (and the previous discussion on the relative distri- 
bution of H I and 70 /um emission) we draw the following 
conclusions: 

- Most of the detected 70 /im emission is coincident 
with Hi column densities of Nhi > 10^^ cm^^ (with a 
peak around 1— 2x 10^^ cm~^). In the case of Holmberg I, 
the peak in the first histogram bin at low surface densities 
is spurious and a result of the lower signal-to-noise in 
these data. 

- At high column densities (> 2.5 x lO^^cm"^), the 
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Fig. 9. — Solid histogram: distribution of the 70 /^m flux density 
above a threshold of 1.8 MJysr~i (~5(t) as a function of Hi col- 
umn density. Dashed histogram: total distribution of H i column 
densities. 

solid (70 fiui) and dashed (H i) histograms follow each 
other closely, implying that most of the high H l column 
density regions are associated with dust emission. In 
other words, there appear to be only few regions of high 
H I column densities that are not associated with warm 
dust emission. 

3.2. IRS spectroscopy 

Selected regions in some of the galaxies in this sample 
have also been observed with the Spitzer IRS spectro- 
graph (Houck et al. 2004) as part of SINGS. In the follow- 
ing we will discuss IRS 5-14/im, low resolution spectra 
obtained for specific regions in the two most active galax- 
ies in our sample, IC 2574 and Holmberg II (Fig. 10, both 
SFR'^O.l Mq yr^^). The spectra have been extracted 
over a circular aperture of a diameter of ~14" that was 
centered on the brightest region seen at 8/im in both 
galaxies (see the caption of Table 2 for coordinates). For 
comparison, Table 3 also summarizes some of the IRAC 
flux densities obtained for the same apertures. 



3.2.1. IRS spectrum of IC 2574 

The spectrum of IC 2574 (upper panel in Fig. 10) 
clearly shows the; broad emission features from polycyclic 
aromatic hydrocarbons (PAHs) which are typically found 
in mid-infrared spectra of massive star-forming galaxies 
(e.g. Telesco 1988, Lu et al. 2003, Smith et al. 2006). To 
derive a simple template for the aromatic/PAH features, 
we have taken the IRS spectrum of the SINGS galaxy 
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Fig. 10. — IRS spectra of IC 2574 (top panel) and Holmberg II (bottom panel). The vertieal dashed line indicates the border of SLl 
and SL2 modules of the IRS spectrograph. The thick grey lines represent the sum of our scaled PAH template (derived from SINGS 
observations of NGC 7552) and the underlying continuum (dashed curve, see text for details) . The black horizontal line indicates the width 
(>10% response) of the IRAC band 4 filter. 



NGC 7552 (12+log(0/H)-8.5, Moustakas ct al. 2007, 
SFR~7M0yr~^) which shows a typical spectrum with 
one of the highest signal to noise ratios in the SINGS 
sample (Smith et al. 2006). From this spectrum we sub- 
tracted a 200 K blackbody spectrum that has been nor- 
malized to the pseudo-continuum at 10 /im and 13.5 /im. 
This PAH template was then scaled to fit the PAH fea- 
ture at 11.3 fim in IC 2574 after a scaled blackbody spec- 
trum of 200 K (shown as the dashed thick line) has been 
added (thick grey line, normalized in the same way as 
for NGC 7552). Although a blackbody is an unphysical 
representation of the dust continuum beneath the emis- 
sion bands, this simple procedure allows us to compare 
the strength of the PAH bands in different sources in a 
consistent way. The relatively low signal-to-noise ratio 
of the spectrum prevents a full spectral decomposition 
(e.g., as done by Smith et al. 2006). 



It is interesting to note that this simple template fits 
the spectrum of IC 2574 quite well to first order. Even 
though the S/N ratio is low, the 7. 7 fiui/ 11. 3 fim PAH ra- 
tio in IC 2574 appears to be lower than in NGC 7552 
- band-variations like that are known for many other 
galaxies (e.g., Drainc & Li 2001, Vcrmeij ct al. 2002, 
Cannon et al. 2006, see detailed discussion in Smith et al. 
2006). Engelbracht et al. (2005) showed that there ap- 
pears to be a metallicity threshold of 12-|-log(0/H)~8.2 
below which PAH emission is not detected in galaxies^^. 
The gas phase metallicity of IC 2574 (Table 1) is slightly 
below their threshold the fact that we do see PAH 

^'^ Note that the oxygen abundances used by Engelbracht ct al. 
(2005) were based on a heterogeneous compilation of measurements 
from the literature, whereas the abundances in our study have been 
derived self consistently, and placed on a common abundance scale 
(see Moustakas et al. 2007 for details) 
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emission in IC 2574 is likely not due to metallicity varia- 
tions within the galaxy (as no evidence of significant local 
metallicity variations has been found in dwarf galaxies; 
e.g., Kobulnicky & Skillman 1996, 1997), but rather due 
to to other factors that influence the aromatic feature 
emission (e.g., radiation field, geometry, elemental com- 
position of the ISM). 

Two bright emission lines are present as well in the 
spectra of IC 2574: [Ne II] emission at 12.8 fj,m and [S IV] 
emission at 10.5 /xm. The latter line is not present in 
the template spectrum of NGC7552 and indicates the 
presence of massive stars - this line is usually found to 
be faint in massive galaxies (e.g., see discussion in Rigby 
& Rieke 2004), but bright in the highly ionized gas of 
blue compact dwarfs (e.g., Madden et al. 2006). 

3.2.2. IRS spectrum of Holmberg II 

The IRS spectrum of Holmberg II (lower panel in 
Fig. 10) is markedly different from that of IC 2574; the 
spectrum is dominated by continuum emission, and only 
very faint PAH features are present. The thick grey line 
again shows the comparison to our simple template: Here 
two blackbody curves (one at 300 K and one at 700 K, 
contributing equal flux density at 8.3 /im) were needed to 
adequately fit the continuum of Holmberg II (note again 
that these curves do not represent physically meaningful 
numbers but are only used to describe the shape of the 
contimmm to first order, see also discussion in Smith et 
al. 2006); we then added a scaled version of our PAH 
template to this continuum emission (thick grey line) . A 
comparison to this curve shows that, although the S/N 
is low, faint PAH features appear to be present at 6.2, 
7.7, 11.3 and 12.7/im- e.g., the 11.3/Ltm PAH feature is 
detected at 7cr (total flux; peak: 4(t). From this it is 
clear that the PAH-to-continuum ratio in Holmberg II 
is much lower than in the case of IC 2574 as discussed 
below. We note that the gas-phase metallicity of Holm- 
berg II (Tab. 1) is lower than in IC 2574 by nearly a fac- 
tor of two (below the threshold derived by Engelbracht 
et al., 2005). This result emphazises the fact that the 
strength of the PAH features is not a simple linear func- 
tion of metallicity (see also Smith et al. 2006, who find 
a wide range of PAH strengths (factor of ~ 10) near the 
Engelbracht et al. threshold). As in the case of IC 2574, 
line emission from [Ne II] and [S IV] is detected in Holm- 
berg II. 

3.2.3. PAH-to-Continuum Ratios 

Using our simple decomposition of the PAH features 
and the continuum emission, we can now constrain the 
PAH-to-continuum ratios for both galaxies. We do this 
for two bands: A) the IRAC band 4 ('8/xm band'), en- 
compassing the broad PAH features at 7.7/xm and 8.6/xm 
(see horizontal line in Fig. 10 for the wavelength range 
covered by the 10% response of the IRAC band 4) and 
B) the PAH band at 11.3/im (here integrated within a 
mock square filter between 10.8 and 11.8/im). We also 
calculate the PAH- to total IR luminosity ratios for both 
regions below (but note that the method employed here 
is not directly comparable to the one used in Smith et 
al. 2006). 

IC 257^: From the IRS spectrum we derive a flux den- 
sity of 2.7 mJy for the IRAC band 4 bandpass which is 



in good agreement with the value derived from the ac- 
tual IRAC band 4 measurement (3.0 mJy, see Table 3). 
The continuum contribution is 0.7 mJy and the con- 
tribution from the PAH feautures is 2.0 mJy; i.e. we 
derive a PAH-to-continuum ratio for this spectral re- 
gion of ~2.9. For the 11 /zm feature we derive a ra- 
tio of 0.5 (continuum: 3.3mJy, PAH: 1.7mJy). Using 
our definitions for the 8/im and ll/zm bandpasses and 
the total infrared (TIR) luminosities in this aperture 
(Ltir, derived from the MIPS images and using the re- 
lation in Dale & Helou 2002) we get the following ratios: 

LpAH,8A<m/LTIR =0.0089, LpAH,llMm/LTIR=0.0011. 

Holmberg II: From the IRS spectrum we derive a flux 
density of 7.9 mJy for the IRAC band 4 bandpass which 
is in good agreement with the value derived from the 
actual IRAC band 4 measurement (8.53 mJy, Table 3). 
The continuum contribution is 7.2 mJy, the contribution 
from PAHs is 0.7 mJy, leading to a PAH to contimmm 
ratio of 0.1. For the 11 /im feature the ratio is even 
lower, 0.02 (continuum: 10.4mJy, PAH: 0.2 mJy). The 
corresponding ratios compared to the total infrared lu- 
minosity in this aperture are: LpAH,8/^m/LTiR=0.0037, 
LpAH,ii;im/LTiR=0.0002. Given the faintness of the PAH 
features in Holmberg II, these values are uncertain (by 
~50%). 

For comparison, we also derive the PAH-to-continuum 

ratios for our template galaxy NGC 7552 and get flux 
density ratios of 6.3 and 0.9 for the IRAC4 and [10.8; 
11.8] /um bandpasses, respectively. For the TIR lu- 
minosities we get the following ratios for NGC 7552: 

LpAH,8A<m/LTIR=0.070, LpAH,llMm/LTIR=0.0077; i.e. 

the ratio in IC 2574 is about a factor of 7 (Holmberg II: 
factor of >20) less than in our template galaxy (see Smith 
et al. 2006 for variations within the SINGS sample). 

3.3. Global relations for the M81 group dwarfs 

In the following we will compare the properties found 
for the M81 dwarf irregular galaxies to other galaxies in 
the 57A/'G5 sample (Kennicutt et al. 2003). To do so, we 
have divided the iS/iVCS' galaxies into four categories: the 
M 81 group dwarf irregular galaxies of this study (shown 
as filled circles in the following plots) , other dwarf galax- 
ies (open circles), elliptical/SO's (open squares) and spi- 
ral galaxies (crosses). All MIPS flux densities for the 
SINGS galaxies are taken from Dale et al. (2006), and 
the gas -phase metallicities are taken from Moustakas et 
al. (2007, based on the Pilyugin & Thuan 2005 strong- 
line abundance calibration). 

As there apparently exists some correlation between 
the Hi and the 70 /xm emission (see discussion in the 
previous sections), we start by comparing the Hi masses 
to the 70 /im specific luminosity^*^ in Figure 11a. It is ob- 
vious that there is a large scatter between the two quan- 
tities (the dashed line indicates a linear relationship). If 
we plot the absolute blue magnitudes (Mb, taken from 
Moustakas et al. 2007) as a function of the 70 /v,m specific 
luminosity instead (Figure lib), this relation gets tighter 
(dashed line). The larger scatter in the Hi 70/im rela- 

In the following wc plot the specific luminosities (units of 
WHz~^). To derive the luminosities in a given MIPS band, this 
number needs to be multiplied by the effective bandwidth of the 
MIPS filter (in Hz). Note that some authors (e.g. Calzctti ctal. 
2005) define the luminosity as the specific luminosity times the 
observed frequency (in Hz). 
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TABLE 2 

IRS vs. IRAC FLUX DENSITIES 



Galaxy 


F(8/im) 


F(6/im) 


F(3.5/im) 


F(8Aim)iRs/F(8/im)iRAC 


F (6/im) IRS /F (6/im)iRAC 




mjy 


mJy 


mJy 






IC 2574=- 


3.10 


1.39 


0.73 


0.89 


0.95 


Holmberg IP 


8.53 


5.75 


2.02 


0.93 


0.81 



^Central coordinates (J2000.0) and diameters D of circular apertures: IC 2574: RA: 10 28 48.3, DEC: +68 28 03, D=13.67", Holmberg II: RA: 
08 19 12.8, DEC: +70 43 08, D=14.63". 
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The dashed line 



H I mass as a function of 70 jum specific luminosity, (b): Blue magnitudes (Mb) as a function of 70 /im specific luminosity, 
in both plots represents a linear relationship between the two variables. 



tion is due to the fac;t that the dwarf irregular galaxies 
of our sample have more Hi mass per blue magnitude 
compared to more massive spirals. This fact, i.e. that 
Mhi/Mb increases for dwarf irregular galaxies was noted 
long before (e.g., Skillman 1996). 

As a next step, we plot the star formation rates (SFRs) 
of individual galaxies as a function of the 70 fim specific 
luminosity (Fig. 12a). The SFRs were derived using SFR 
= L(Hacorr)/l-26xlO'*icrgs-i [Mq yr^i] (Kcnnicutt 
1998) (assuming solar mctallicity and the Salpctcr (1955) 
initial mass function between 0.1 and lOOM©). L(HQ;corr) 
is the Ha luminosity corrected for extinction within the 
galaxy. Here we use the relation derived for SINGS galax- 
ies by Kcnnicutt ct al. (2006) and Calzctti ct al. (2006): 
F(Ha)corr=F(Ha)obs+0.035xF'(24/im). F(Ha)obs is the 
observed Ha flux (in ergs~^cm~^, corrected for Galac- 
tic absorption and contribution from [N H]) taken from 
Kcnnicutt et al. 2006 and Lee 2006; F'(24/im) is defined 
as F'(24/xm)=1.25xl0i3Hzx F(24/xm)[Jy] (i.e. F(24y^m) 
multiplied by the observed frequency in Hz). SFRs were 
only derived for those SINGS galaxies for which accu- 
rate Ha measurements are available. The tighter relation 
between L(70/im) and SFR (as compared to the Mb- 
L(70/xm) relation discussed above), can (to first oder) be 
explained by the fact that, typically, 70/xm is close to the 
peak of the FIR emission and thus is a measure of the to- 
tal bolometric luminosity which is expected to scale with 
the SFR. 

However, from Figure 12(a) it is also obvious that the 



galaxies do not follow a linear relation over all luminosi- 
ties (as indicated by the dashed line). To exemplify this, 
the ratio of the 70 /xm specific luminosity and the SFR 
is plotted in Figure 12(b) as a function of the gas- phase 
oxygen abundance (taken from Moustakas et al. 2007). 
If there was a linear relation between the 70 /im emis- 
sion and the SFR, one would expect the symbols to lie 
near a horizontal line in this plot. However, the dwarfs, 
shown as circles, clearly occupy a lower Lyo^^m/SFR space 
compared to the more massive galaxies. On average, 
the dwarf galaxies appear to be underluminous in 70 fim 
emission relative to their SFR by a factor of ~ 2. A sim- 
ilar conclusion (dwarfs have lower Lir/Lho) has been 
reached by Hunter et al. 1989, albeit for more luminous 
systems. This behavior may, to first order, be attributed 
to the fact that these objects have low metallicities and 
low dust contents (see discussion in Sec. 3.5 and Draine 
et al. 2007). However it is also clear that the scatter is 
large and that there is no simple relation between oxygen 
abundance and 70/im luminosities per unit SF. 

We now investigate how the global MIPS colors (here: 
the 70/um/160 /zm and 70/um/24/um ratios) of the M81 
groups dwarfs compare to the other SINGS galaxies. 
In Figure 13 we plot both ratios as a function of the 
oxygen abundance: Although the scatter is large, the 
dwarf galaxies have elevated 70/im/160/im ratios and 
70/zm/24 /um ratios compared to the spiral galaxies in 
SINGS. The elevated 70/im/160/xm ratios imply that 
the effective temperature of the dust in the dwarf irreg- 
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ular galaxies is on average higher than in more massive 
spirals. In this simphstie pieture, the higher effective 
temperature results in a peak of the SED that is shifted 
toward the 70 /xm waveband (cf. Hunter etal. 1989, Dale 
ct al. 2005). In this context it is interesting to keep in 
mind that the L70prn/SFR ratio in dwarfs is lower than 
in the spirals (see above). In other words, if the dust 
temperature in our sample dwarfs were the same as in 
the spirals, the Lyo^m/SFR ratio would decrease even 
further. 

The elevated 70 /j.m/24 jim ratio in our sample dwarfs is 
more difficult to interpret as the origin of the 24/im emis- 
sion is not certain. According to the models by Draine 
& Li (2006), this emission is due to in part to single- 
photon heating, although in galaxies with strong 24/im 
emission it is primarily due to warm grains in strong ra- 
diation fields. The 24/im luminosity thus depends on the 
intensity of the radiation heating the dust (which de- 
pends on the density in the HII regions, the degree of 
clustering of O stars, as well as on the dust abundance) . 
Future detailed modelling of the SEDs is needed to fully 
describe this behavior. We also note that our results 
do not necessarily hold for all classes of dwarf galaxies; 
i.e., extreme cases such as the metal-poor blue compact 
dwarf SBS 0335-052 have more extreme colors (very low 
70/tm/24/tm ratio, Houck et al. 2004). 

3.4. Spatial comparison to Ha 

The apparent strong correlation between the 70 /im lu- 
minosities and the SFR implies that the ongoing star 
formation is the main heating source for the warm dust. 
We now compare the spatial distribution of the Ha and 
the 70 /im emission: Figure 14 (top) shows the Ha im- 
ages for the three brightest Ha emitters in our sample 
(IC 2574, Holmberg II and Holmberg I). The contour 
shown in the invidual panels represents low-level 70 /tm 
emission (shown in Figure 14, bottom). From this we find 
a very good correlation between the locations of strong 
Ha and 70 /tm emission. This has been found in other 
galaxies before (e.g., M51: Calzetti etal. 2005) but it 



is interesting to note that the same holds true even for 
faint dwarfs such as Holmberg I. We note that there is 
also diffuse dust emission present in our objects (e.g., 
Holmberg II, IC 2574) that is likely heated by the un- 
derlying stellar population or by UV photons that are 
leaking from the HII regions (see Cannon et al. 2006b for 
a detailed discussion on the diffuse dust component in the 
Local Group galaxy NGC 6822, see also, e.g., Popescu et 
al. 2002, Popescu & Tuffs 2003, Hinz et al. 2006). 

3.5. Dust Masses 

We now discuss the dust masses of the M81 group 
dwarf irregular galaxies. As the dust mass derivations 
mostly rely on the longest wavelength MIPS bands (i.e. 
the 70 /im and 160 /tm bands), and given the large error 
bars in the measurements of the dwarf galaxies discussed 
here (in particular in the 160 /im band, see Table 1), the 
dust masses presented in the following can only be con- 
sidered to be order of magnitude estimates. Also, the 
observations are not sensitive to cold dust, which emits 
at wavelengths >>160/im (i.e., there may be more dust 
present that is not heated by star formation). 

To estimate the dust masses we use the new model pre- 
sented in Draine & Li that include variable PAH abun- 
dances (2006, which is based on the models presented 
in Li & Draine (2001, 2002)). In their model, radia- 
tion field strengths are varied via power-law distribu- 
tions; PAH, silicate and graphite grains are illuminated 
and the resulting SEDs can be compared to the observa- 
tions. Table 3 (column 2) summarizes the dust masses of 
the sample dwarf galaxies using their models (note that 
DD0165 and M81dwA were non-detections; the flux 
density measurements of MSldwB are too uncertain to 
derive a meaningful dust mass): the dust masses span 
almost two orders of magnitude from '^l.O 70x10^ Mq. 
If we use the SED models by Dale etal. (2001) and Dale 
& Helou (2002) we get masses that are higher by a fac- 
tor of '--^5. We attribute the difference to the fact that 
their models and dust mass correction factors were de- 
veloped for normal galaxies and are therefore likely not 
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TABLE 3 
Dust and Gas Masses 



Galaxy'' 


(10* Mq) 


Mhi 
(10« Mo) 


(lO^Mo) 


Mdust/Mni 
10-3 


10-3 


IC 2574 


72 


14.75 


2.8 


0.49 


2.6 


Holmbcrg II 


12 


5.95 


1.2 


0.20 


1.0 


Holmberg I 


6.8 


1.40 


0.4 


0.48 


1.7 


DDO053 


1.0 


0.60 


0.2 


0.17 


0.5 



^See Table 1 for galaxy parameters 

'^Dust mass derived using the SED models of Draine & Li 2006 
^HI mass in region where dust emission is present 
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Fig. 13. — (a): 70 /im to 24 fj,m flux density ratio as a function of oxygen abundance, (b): 70 /um to 160 /im flux density ratio as a function 
of oxygen abundance. 




Fig. 14. — Top: H« images for the three brightest dwarf galaxies in our sample (a: IC 2574; b: Holmberg II, c: Holmberg I). Contours 
indicate 70 fun surface brightness levels [(a): 1.8MJysr-^, (b); 4MJysr-^, (c); 2.6MJysr-^. Bottom: 7Qf^ra maps of the three galaxies 
with the same contour overlaid (the emission seen toward the north-west of Holmberg I is noise) . 
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appropriate for dwarfs. 

With these dust masses wc can now estimate the dust- 
to-gas ratios of the dwarf galaxies in our sample. In the 
following we define Mgas=MHi for the dwarfs, i.e. we do 
not take the possible (unknown) contribution of molecu- 
lar gas into account (and we do not correct for the con- 
tribution of Helium) . Using our dust and H I masses 
(Table 3, column 3), we derive an average Mdust/Mni 
ratio of ~ 3 X 10"'* (numbers for the individual galax- 
ies are given in Table 3, column 4). These values do 
not appear to be a function of the measured metallici- 
ties (however the spread in metallicies in the remaining 
sample is small) . If we consider only the H I mass of the 
M 81 group dwarfs in an aperture defined by the extent 
of the dust emission (M^']*, Table 3, column 5) this ra- 
tio increases to an average of Mjust/M^'^* ~ 1.5 x 10~^ 
(Table 3, column 6). This is because, on average, about 
three quarters of the extended H i mass has no associated 
bright dust emission (cf. Figs. 1-7). 

Draine etal. (2007, in prep.) present detailed dust 
mass estimates for all SINGS galaxies and derive a typ- 
ical spread in dust-to-gas ratios for the more massive 
galaxies of 0.005<M(iust/Mgas <0.02 (here including the 
contribution of molecular hydrogen). The value for the 
dwarfs is thus about an order of magnitude lower; it 
would decrease further if the dwarfs had a significant 
component of molecular hydrogen. The ratio for the 
dwarfs would increase, however, if they had an addi- 
tional cold dust component (e.g., in the outskirts where 
the stellar radiation field is low, see also Draine et al. 
2007) that can not be traced with the available MIPS 
observations. E.g., in the case of NGC 1569, Galliano et 
al. 2003 report a millimetre excess in the dust SED - 
this cold dust component could account for up to 70% 
of the total dust mass in this object. Such a cold reser- 
voir may thus increase the total dust content by a factor 
of a few, but likely not by an order of magnitude. We 
thus conclude that the dust-to-gas ratio (Mdust/Mgos) 
in our sample dwarfs appears to be significantly lower 
than what is found in spiral galaxies (which can not sim- 
ply be explained by a linear scaling of this ratio with 
metallicity, see also the discussion in Draine et al. 2007). 

4. SUMMARY 

We present observations of warm dust and atomic gas 
in seven dwarf irregular galaxies in the M81 group us- 
ing data from both SINGS and THINGS. Five of the 
seven targets have been detected with Spitzer out to 
160 /im (the dwarfs with the lowest star formation rates, 
M 81 dwA and DDO 165, are non-detections). As molec- 
ular gas in these systems has yet to be detected (only 
likely exception: IC 2574, Leroy et al. 2005) the Spitzer 
observations give a first glimpse of the nature of the 
non-atomic ISM in these galaxies and provide impor- 
tant information to design follow-up observations of the 
molecular gas phase (e.g., to select regions of interest for 
pointed observations using millimeter telescopes). 

We find that the warm dust emission as traced via the 
70 /im observations is associated with high H i column 
densities (Nhi ~ 1 x lO^^cm"^), close to the 'canon- 
ical' star formation threshold found by previous stud- 
ies (e.g., Skillman 1996, Walter & Brinks 1999, Schaye 
2004) . Most regions with H i column densities of Nhi > 
2.5 X 10^^ cm~^ have dust emission associated with them. 



For the brightest regions at 70 /im there is a good correla- 
tion with the location of HII regions, indicating that ac- 
tive star formation is needed to heat up the dust locally. 
However, in some cases, there is diffuse dust emission 
present at larger radii which does not appear to coincide 
with compact HII regions. This diffuse emission is likely 
due to the re-processing of non-ionizing photons in the 
ISM or the escape of radiation from the star formation 
regions (see also Cannon et al. 2006b). 

IRS spectroscopy in the brightest regions in IC 2574 
and Holmberg II (which have comparable star forma- 
tion rates) reveal distinctly different spectral shapes: 
Whereas PAH features are clearly detected in the spec- 
trum of IC 2574 those features are weaker in Holmberg II 
(which has lower metallicity) by an order of magnitude. 
This emphazises that the strength of PAH features is not 
a simple linear function of metallicity (sec also Smith ct 
al. 2006). The PAH to TIR contimmm ratio in IC 2574 
is a factor of ^7 less than what is found in a typical 
SINGS spiral (a factor of > 20 in the case of Holm- 
berg II). 

While the H i masses arc well-constrained, it is diffi- 
cult to constrain the dust masses with a high degree of 
certainty. We estimate dust masses of ~10'*-10^ Mq for 
individual targets, resulting in an average dust-to-gas 
ratio (Mdust/Mm) of ~ 3 x 10"^ (1.5 x 10"^ ff only the 
H I that is associated with dust emission is considered) . 
This can be compared to the range in values derived for 
the SINGS galaxies by Draine etal. (2007, in prep.) of 
0.005<Mdust/Mgas <0.02. Thus, the average dust-to- 
gas ratio in the dwarfs is lower by about an order of 
magnitude as compared to more massive spirals (a find- 
ing that can not simply be explained by a linear scaling 
with metallicity, see also Draine et al. 2007). Future sen- 
sitive observations at longer (sub-mm) wavelengths are 
critical to constrain the possible presence of a colder dust 
component in the galaxies (not heated by the stellar pop- 
ulation) that may be missed by the MIPS observations. 

We also find that the dwarf galaxies in our sample are 
underluminous at 70 fj.m for a given SFR by about a fac- 
tor of ~ 2 compared to the more massive and metal-rich 
galaxies in SINGS. However, interestingly, the average 
70 //m/160 /im ratio in the dwarfs is higher (factor of ^2) 
than in the spiral galaxies. In a simplistic picture, this 
can be attributed to higher effective dust temperatures 
in the dwarf galaxies (which shifts the peak of the warm 
dust SED toward 70mm, cf. Dale etal. 2005). Similar 
conclusions on the dust temperature have been derived 
by other authors studying more luminous dwarf systems 
(Hunter etal. 1989, Dale etal. 2005, Engelbracht etal. 
2005, Cannon etal. 2005, 2006a, 2006b) and have been 
explained in the context of stronger radiation fields in 
the dwarfs. It is interesting to note that, if the dwarf 
galaxies had the same temperature as the more massive 
spirals, the 70/im luminosity for a given SFR would de- 
crease further (relative to the spirals). Overall, there is 
a better correlation between the SFR (or optical mag- 
nitudes) and the 70 /zm luminosity than between the H I 
mass and L(70/im). This provides additional evidence 
that the FIR emission in the sample dwarf galaxies is 
powered by ongoing star formation and does not strongly 
depend on the total Hi mass of the galaxy host. 
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